Effects of thermal fluctuation and the receptor-receptor 
interaction in bacterial chemotactic signalling and adaptation 



Yu Shi * 

Cavendish Laboratory, University of Cambridge, Cambridge CBS OHE, United Kingdom 

Abstract 

Bacterial chemotaxis is controlled by the conformational changes of the re- 
ceptors, in response to the change of the ambient chemical concentration. In 
a statistical mechanical approach, the signalling due to the conformational 
changes is a thermodynamic average quantity, dependent on the temperature 
and the total energy of the system, including both ligand-receptor interaction 
and receptor-receptor interaction. This physical theory suggests to biology a 
new understanding of cooperation in ligand binding and receptor signalling 
problems. How much experimental support of this approach can be obtained 
from the currently available data? What are the parameter values? What 
is the practical information for experiments? Here we make comparisons 
between the theory and recent experimental results. Although currently com- 
parisons can only be semi-quantitative or qualitative, consistency is clearly 
shown. The theory also helps to sort a variety of data. 
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I. INTRODUCTION 



Bacterial chemotaxis refers to the phenomenon that a bacterium such as Escherichia coli 
swims towards higher concentration of attractant or lower concentration of repellent |l]-f§]. 
With the switching rate determined by the change of the ambient chemical concentration, 
the motors of the bacterium switch between counterclockwise and clockwise rotations, con- 
sequently the bacterium switches between tumbling and running. The ratio between the 
frequencies of the two rotation modes is determined by the rate at which kinase CheA phos- 
phorylates CheY, which binds the base of a motor. CheA phosphorylation rate is regulated 
by the receptor conformational state, which is influenced by ligand binding. The receptors 
are dimeric and joined to a CheA dimer by a CheW dimer, furnishing a signalling complex. 
Hence a receptor dimer can be regarded as a basic unit, as supported by the finding that a 
receptor dimer with a damaged subunit can still work || . Because of thermal fluctuation, 
even in the absence of ligand binding, or in a fully adapted situation, there is still a certain 
probability distribution for a receptor dimer to be in different conformational states; micro- 
scopically a receptor dimer stochastically flips between the two states. Attractant binding 
changes the probability distribution, causing the receptor dimer to be more likely in the 
state corresponding to the lower CheA phosphorylation rate. On a longer time scale, after 
an initial response to ligand concentration change, the activity of the system returns to the 
pre-stimulus level. A careful consideration of such a basic picture already finds the ideas of 
statistical mechanics necessary: with the presence of thermal fluctuation, it is the proba- 
bility distribution of the conformational states of the receptors that is monitored by ligand 
concentration change and determines the motor rotation bias. However, it seems that this 
point is not universally appreciated in biological literature. 

The chemotactic response is very sensitive ||. It had been conjectured that there might 
be cooperation between receptors or the signalling complexes so that the signal could be 
amplified |7],[3|]. The fact that most of the receptors cluster together at a pole of the cell 
provides further clues for cooperation between receptors More importantly, it was 



2 



found experimentally that the clustering of receptors is not favorable for counting statistics 
and that the receptor cluster does not favor a special end of the cell ||10|| . This is an indication 
that there is a special reason for clustering, which may well be to make the receptor-receptor 
interaction possible. 

With a detailed analysis on the experimental findings, I suggested the possible existence 
of interaction between neighboring receptor dimers and constructed a statistical mechanical 
theory to provide a picture of how the receptors cooperate through physical interaction 
and how the thermal fluctuation makes statistical mechanics important in the signalling 
process [11,|12[. In our model, we combine cooperativity and feedback to account for the 



sensitivity and adaptation. As will be stressed here, the first message from our approach is an 
emphasis on thermal fluctuation. In a cell, the energy scale is comparable with the thermal 
fluctuation. Moreover, thermal fluctuation helps to distinguish different stimuli. Because of 
large separation of time scales, the thermal fluctuation can be treated as quasi-equilibrium, 
so equilibrium statistical mechanical can give a reasonable response-stimulus relation. Hence 
the basic elements of our theory is useful no matter whether there is interaction between 
receptor dimers. The second message of our theory is that the anticipated cooperation is 
just physical receptor-receptor interaction between neighboring receptor dimers. Therefore 
the conformational state of a receptor dimer is not only influenced by ligand binding of 
itself, but also by the receptor-receptor interaction which depends on the conformations of 
the two neighboring receptor dimers. The third message is that the large separation of time 
scales leads to a complementary usage of equilibrium statistical mechanics in calculating the 
response on a shorter time scale and a non-equilibrium description of the adaptation on a 
longer time scale. Dynamics on the longer time scale determines whether randomness of 
ligand binding is quenched or annealed on the shorter time scale of quasi-equilibrium state, 
as will be elaborated later on. In the high temperature limit, this does not make a difference 
on the average signalling. 

Recently there appeared some experimental data which are more directly relevant for the 
many-body nature of the receptor cluster and the possible cooperation []T3| P~5p . Therefore it 



is interesting and important to make comparisons between the theory and the experimental 
results, testing the theory on one hand, and providing some information on what experiments 
are wanted on the other hand. However, we do not expect the model in current form to 
perfectly fit everything on this complex system. Rather, what we provide is a theoretical 
framework, on which refinements are possible. For example, we have only considered the 
cooperation between the receptor dimers, while extensions to possible cooperation among 
other components at later stages of the signalling process, for example, CheA, CheY, CheZ 
and the switch complex, is straightforward if sufficient information is available. The idea 
of receptor-receptor interaction broadens the view on cooperation, which previously largely 
meant the existence of more than one binding sites, as described by the model presented by 
Hill a century ago ||16|| . For simplicity, we try to preserve the scenario of one binding site, 
while extension to the situation with more binding sites is straightforward if necessary. Our 
strategy is to start with the minimum model, which yet explains the most essential features. 

With improvement and simplification, we first describe the theory. Then we make com- 
parisons with the experimental results, followed by summary and discussions. 

II. THEORY 

Consider a lattice of receptor dimers, as shown in Fig. 1. Let the coordinate number 
be z/, which is 6 for a honeycomb lattice and is 4 for a square lattice. The exact coordinate 
number in reality is subject to experimental investigations. The behavior of the system is 
determined by its energy function, or Hamiltonian, which can be written as 

Hit) = - £ TijVtVj - £ HiVi + £ WiVi. (1) 

<ij> i i 

Vi is a variable characterizing the conformation of receptor dimer i. It may be interpreted as 
the position of the receptor molecule with respect to a certain equilibrium position. In the 
popular two-state approach, Vi assumes one of two values V° or V 1 . Hi is the influence, or 
force, due to ligand binding and the modulation of methylation level. Hi = if there is no 
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ligand binding, while Hi = H if there is a ligand binding. —HiVi is the energy due to ligand 
binding, hence ligand binding causes the energy difference between the two conformations to 
make a shift of HiV 1 — V°). Wi(V° — V 1 ) is the original energy difference between the two 
conformations, (ij) denotes nearest neighbouring pairs, —TijViVj is the interaction energy 
between the neighboring receptor dimers. 

For convenience, defining Si = 2(Vi — V°)/AV — 1, with AV = V 1 — V°, one transforms 
the Hamiltonian to 



where Si = 1,-1 represents the two conformational states of the receptor dimer at site i, 
Jij = T tj AV 2 /A, Bi = HiAV/2, U t = AVWi/2- AV 2 EjTy. We refer to B { as "field". For 
simplicity, it is assumed that = J and U{ = U are independent of % and j. Bi = if there 
is no ligand binding, while Bi = B = HAV/2 if there is a ligand binding. Hence energy 
difference due to ligand binding between the two conformations are 2B. USi represents the 
original energy in the absence of ligand binding. Eq. ([!]) and (0) can be justified as follows. 
It is reasonable to assume an interaction energy proportional to (Vi — V,) 2 , which can be 
reduced to —T^ViVj, with constant terms neglected and the terms proportional to Si or Sj 
included in J^iUiSi. This assumption is simple enough to allow a feasible treatment which 
yet captures the essential features. 

From now on, we focus on Eq. (13). Suppose that before time t = 0, there is no ligand 
binding the system, or the system is fully adapted though it is bound to ligands. Hence 
Bi(t < 0) = 0. Afterwards, at time t — 0, the occupancy, i.e. the fraction of receptor dimers 
with ligands bound, changes to c. Hence the occupancy change is 5c = c. This means 



0, with probability 1 — c 
The occupancy c is determined by the ligand concentration L, c = L/(L + K4), where the 
dissociation constant is on a time scale during which the receptor has undergone many 



Hit) 



^2 JijSiSj — ^2 Bi(t)Si + ^2 UiSi 



(2) 



Bi{t = 0) = BP, with 
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flips between different conformations, hence it is an average and phenomenological quantity. 

On the other hand, through the modulation of methylation level by CheB and CheR, 
there is a negative feedback from the receptor state Si to the field Bi, with a time delay t r . 
A simple quantitative representation of this feedback is 

dBi(t) 



dt 



-a[Si(t-t r ) -mo], (4) 



where a > 0, m is the pre-stimulus average of Si. Precise forms of both the energy function 
and the feedback are, of course, subject to experimental investigations. It seems that in 
biological world, feedback is a ubiquitous way to achieve adaptation and preserve sensitivity 
of response. 

A remarkable feature of this system is the large separation of time scales. Ligand bind- 
ing and conformational change occur within only millisecond, while overall time needed to 
complete the adaptation, through the slow modulation of methylation level, is on the time 
scale of many seconds to minutes [|T8|,3 . We note that in most cases, ligand debinding is on 



a much longer time scale than ligand binding, seen as follows. Consider the kinetics of the 
following reaction 

L + R^R L , (5) 

where R represents the receptor without ligand binding, while Rl represents the liganded 
receptor. k + and /c_ are reaction rates for the binding and debinding, respectively. The 
ratio between the time scales of debinding and binding is k + L/k^ = L/K^, where is the 
dissociation constant. A typical value is Kd ~ 1.2/iM [[|. Usually, L is much larger, so 
the debinding time scale is much longer than the time scale of ligand binding and receptor 
conformational change. In extreme cases, when L is comparable with Kj, debinding time 
scale is comparable to binding time scale. 

With the large separation of time scales, the treatment within the above framework 
becomes easier. One may discretize the time on the scale of adaptation, according to the 
feedback delay time, t is thus replaced by an integer r, which is the integer part of t/t r . On 



the other hand, each instant r is still very long compared with the time scale of conformational 
change. Hence the activity at each r is an average quantity m(r), which can be calculated 
from the Hamiltonian in (|2|) by standard methods of equilibrium statistical mechanics. The 
average activity m is just on the time scale of the measurement of the macroscopic quantities 
such as motor bias, longer than the very short period in which the receptor is in either of 
the two conformations, but shorter than the adaptation time. In making the average, an 
important thing is that the randomness of the field is usually quenched since L » Kd, and 
is annealed otherwise. In fact we obtain a generalized version of the so-called random-field 
Ising model; in the conventional random-field Ising model, the average field vanishes, but it 
is generically non-zero in our model. On the long time scale, the field changes because of 
feedback. It can be expressed as B^t) = B® + M(r), where M(r) is an induced field due 
to methylation modulation, 

M(r) = -aJ2[m(k)-m ]. (6) 

fc=0 

Before stimulation, m(r < 0) = mo is determined by U . If and only if U — 0, mo = 0, 
which means that each receptor dimer is in either of the two conformational states with 
equal probability. 

In most cases, the randomness of Bf is quenched, the general relation between m(r) and 
Sc is then 



m(r) 



28c 



l+cxp[-2/3(i/Jm(T)-0(r-l)o- V (m(k)-m )+U+B)} 

j 2(l-fc) 1 / 7 n 

l+cxp[-2/3(i/7m(r)-0(r-l)o- V [_ 1 (m(k)-m )+U)} ' V > 



where j3 = 1/ksT, the step function 8(x) is 1 if x > 0, and is otherwise. On the other 

hand, when ligand concentration is comparable with Kd, the randomness of Bf is annealed. 

Then it can be found that 

§ c [ e p(f(m)+B) _ e -/3(/(m)+s)] + (i _ 8c)[ePK m ) - e^" 1 )] 
m = ^pC/M+s) _|_ e -/3(/(m)+B)] + (1 - 5c)[e^/( m ) + e"^/( m )] ' ^ 

where f(m) = vJm — 6{r — l)cY^k=o m W + ^ ■ 



m(r) vs. c relation corresponds to the response-stimulus relation. After the step increase 
at r = 0, m(r) always decreases towards the pre-stimulus value m®. This explains the 
robustness of exact adaptation ||19|| . In practice the adaptation time is obtained when m-mo 



reaches the detection threshold m*. 

The results can be simplified under the condition that the thermal fluctuation is so strong 
that f3vJ and j3B are not large. Then both Eq. (0) and Eq. (||) can be simplified to 

(3B5c f (3a V 

with 

m ° - T*h- (10) 

1 — (3vJ represents the enhancement of response compared with non- interacting scenario. 

One may obtain the adaptation time t*, after which m — mo is less than the detection 
threshold m*: 



t _ log 5c + log(j5j) - log m* 



m* can be related to the lower bound of detectable occupancy change, 5c*, by 



(11) 



m = T=fcj> (12) 



hence 



r* = lQgfc - lQg / C \ (13) 

-Mi -4b) 1 ) 

At exact adaptation, setting m(r) = mo, one may obtain the total induced field due 
to methylation modulation as M* = Be. Then for the next stimulus, suppose that the 
occupancy changes from 5c to 5c + Ac at a later time r 1; it can be found that the result 
with the occupancy 5c + Ac and the induced field M* is the same as that for the situation 
in which the occupancy is Ac and there is no induced field. That is to say, the previous 
occupancy change has been canceled by the induced filed M*, therefore the fully adaptation 
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with ligand binding is equivalent to no ligand binding. So m(r > t%) is given by the above 
relevant equations with r changed to r — r 1; and 5c substituted by Ac. One can thus simply 
forget the pre-adaptation history, and re-start the application of the above formulation with 
T\ shifted to 0. The cancellation holds exactly only under the assumption of small j3uJ 
and f3B, which is likely the reality. The finiteness of detection threshold further widens the 
practical range of its validity. 

III. COMPARISONS BETWEEN EXPERIMENTS AND THE THEORY 

A. Clustering. 



The clustering has recently been experimentally studied in greater detail [15]. The ob- 
served clustering of receptors and the co-localization of the CheA, CheY, and CheZ with 
the receptors is a favor for the effects of interactions. An in vitro receptor lattice formation 
was also observed |T7| ]. 

B. Response-stimulus relation. 

A basic prediction of our theory is the response-stimulus relation. Note that the time 
scale of the response, corresponding to m in our theory, is longer than the very short life- 
time of a specific conformation, but is only transient on the time scale of the adaptation 
process. An interesting thing is that m in our theory is measurable. Motor rotation bias 
was measured |13| . From this result we can obtain m, as follows. The motor bias is 



b fccw/(fccw fcw)i (14) 

where f ccw and f cw are rates of counterclockwise and clockwise rotations, respectively. Sup- 
pose the value of b is r% for conformational state 1, and is r_i for conformational state —1. 
Then the average bias is: 

b — rix + r_i(l — x), (15) 



where x is the average fraction of receptors with state 1. x is related to m by m = x—(l—x) = 
2x — 1. So if we know r± and r_i, we can obtain m from b. However, there seems to be 
no investigation on r% and r_j. A simple assumption which is often implicitly assumed is 
that ri = 1, r_i = 0, i.e. state 1 corresponds to counterclockwise rotation, while state —1 
corresponds to clockwise rotation. We follow this assumption here. But it should be kept in 
mind that an experimental investigation on r% and r_i would be very valuable. Therefore, 
for the time being, we use 

b = ^—> ( 16 ) 

Thus from the pre-stimulus value of b, one may determine mo, and thus (3U. An empirical 
formula is b = 1 — 0.0012(rcc? — 360), where red is the absolute angular rate of change of 



direction of the cell centroid in degree ■ s [13,23]. From [23|, the pre-stimulus value of red 



is known as ~ 600, so the pre-stimulus value of b is ~ 0.712. Hence 

mo = z ~~n 7 ~ 0.424. (17) 
1 — pz/J 

The occupancy change used in [T3|] was calculated from the ligand concentration under 
the assumption that the ligand randomly binds one of two possible binding sites: in addition 
to the site with Kd ~ 1.2/iM, as widely acknowledged ||18|| , there is another site with 

~ 70/iM. This was based on an earlier attempt to have a better fitting for the adaptation 
time P0|| . However, as said above, we try to make things as simple as possible in the first 
instance, so prefer to preserve the scenario of one binding site with ~ 1.2/zM. Actually 
with one binding site, as discussed later on, it seems that our theory can fit the adaptation 
time by choosing appropriate parameter values, thus improve the coherence among various 



data. So we should first transform the values of the occupancy change in [13| to the values 



one would have obtained without the assumption of two binding sites. One has 

cj = ^(c 1 + c 2 ), (18) 

where cj represents the occupancy used by Jasuja et ai, C\ corresponds to dissociation 
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constant K\ = 1.2/iM, c 2 corresponds to dissociation constant K 2 = 70/iM. From q = 
Lj (L + Ki), 1 = 1,2, one obtains the change of the occupancy 

K16L 

° Cl (L + 5L + K{) (L + Ki) ' U Jj 

where <5L is the change of ligand concentration. Since 5L << L, one may obtain 5c\ = 
25cj/(l + a), where a w + K^ 2 / K 2 (L + fsT 2 ) 2 . With L w 10/iM, awl, one has 



5c\ w 5c j. Therefore under this condition, we may simply use the occupancy change in |L3 
Eq. (|T6|) leads to the relation between the initial change of m and that of the motor bias, 
5b, 

5m = 25b, (20) 

where 5m = m(5c, r = 0) — m . 

So the data in Fig. 3 of [I3~| can be transformed to 5m vs. 5c relation as shown in our 



Fig. 2. Unfortunately, it is notable that the data is limited to very low values of occupancy 
change! Nevertheless, a qualitative fitting can be made. According Eq. (H), setting r = 0, 
we fit the data with a straight line 5m = a5c, where 

(3B 



1-puJ 



10.49 (21) 



is the slope of the fitting line. 



C. Adaptation time. 

Eq. (|I~9"D tells us that with a same concentration change, the occupancy change and 
thus the response decreases with the increase of pre-stimulus ligand concentration. This 
is verified by Fig. 7 of ||20|| . Eq. (|TT| ) predicts that the adaptation time increases linearly 
with, but not proportional to, the logarithm of occupancy change. It had been thought that 
the adaptation time is proportional to the occupancy change [^T|,^,^]. We found that a 
logarithmic relation is also consistent with the currently available data. As an example, we 



examine the better set of the data, the left plot (D-ribose), in Fig. 4 of p2 |. For accuracy, 
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the two data points at the highest and lowest concentration changes are dropped. This 
is because they are at the detection limits, and they have no recognizable differences in 
adaptation time with the data points closest to them respectively, although the values of 
concentration change are quite different. Moreover, the adaptation time is recorded to be 
zero for the two smallest values of concentration change, so the data point with the smallest 
concentration change should be ignored. Using Kj, = 3 x 10~ 7 (no unit was given, but 
should be the same as that of the concentration, so there is no problem in using it), we 
transform the concentration to the occupancy. The transformed data is shown in our Fig. 3, 
with Fig. 3a the normal- normal plot and Fig. 3b the normal- log plot. While there could 
be a fitting with a proportional relation, as usually assumed, it is at least reasonable to 
fit them with a logarithmic relation, t*/0.354s = r* • t r /0.354s = (g log 10 5c + /i)/0.354s, 
with g = 95.151 x 0.354s = 33.7s and h = 124.0574 x 0.354s = 43.9s. The factor 0.354s 



comes from the data normalization in |p2]| , which is the percentage of one of three maximum 
recovery times, 0.56m, 0.58m and 0.62m, i.e. 35.4s on average. From Eq. (|Tl"|), we have 



-^(1-4; 

and 

t r [log 10 6c*} 



Pa ,-9- (22) 



(23) 



Using 5c* ~ 0.004 p0| , and assuming t r ~ 0.1s, one finds 

— ^ — « 0.0068 to 0.013. (24) 

where the first value is estimated by using (0), and the second by using (|23|). They are 
close to each other, as an indication of the consistency of the theory. 

Furthermore, our predicted logarithmic relation may explain the discrepancy in the anal- 
ysis of the data in Fig. 4 of [|20| about a relation between the adaptation time and the 
concentration. The logarithm can simply decrease the predicted value of adaptation time, 
without resorting to the assumption of the existence of two binding sites. We have tried to 
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make a quantitative fitting for the data in Fig. 4 of [ 20| . Using Kd = 1.2//M, we transform 



the ligand concentration to the occupancy, as shown in our Fig. 4. To make better use of the 
data, we ignore data point for 5c > 0.95, because the finiteness of detection threshold may 
cause uncertainty in deciding the adaptation time; the data for 5c > 0.95 show too large 
variation for so close values of 5c. The fitting straight line is t* = r* • t r = g log 10 5c + h, 
with g = 156.3513 and h = 114.9912. Using (|22]) and ([23|), and same values of 5c* and t r as 



above, one finds 



- ■ 0.0015 to 0.0047. (25) 



Again, these two numbers are close to each other. Moreover, ( ^4|) and ( ^5|) are of the 
same order of magnitude, though they are obtained from different sets of data. Even closer 
numbers may be obtained by tuning the vale of t r . 

D. Che A activity. 

Another interesting and important experimental result is on the relative CheA activity, 



which has been analysed by using Hill model with a non-integer coefficient ||14|| . Here we 
examine the data from the perspective of our theory. 

Suppose S = 1,-1 correspond respectively to CheA activity Ai and A-\. Then the 
average CheA activity is \{A\ + A_i) + ^{A\ — A_i). Consequently the relative CheA 
activity, as measured in fnfl , is 

= (A. + A^ + ^-A^mjSc) _ L 

(A 1 + A. 1 ) + (A 1 -A^ 1 )m(5c = 0) L + K d ' 1 ' 

where F = -^V, with E = + Al/(A- 1 - Al) > 0. Note that A- X > A- It 



B 



is constrained that for attractant binding, F < 1, since R > 0. Setting F = 0.95 and 



Kd = 20/iM, we obtain a reasonable fitting to Fig. 1 of fl4| , as shown in our Fig. 4. 
Therefore 

E ^a(—--), (27) 

v 0.95 B !l v ; 
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which, combined with Eqs. (|i~7|) and (|21|), implies that the ratio between the two levels of 
CheA activity is A^xj A\ « 164.77. Very interestingly, this result of deduction is in good 
consistency with the available experimental information that this ratio is more than 100 [[|. 
Again, this is an indication of the consistency of the theory. 

We note that there is discrepancy in the fitting. This may be because of some other 
factors not considered here, especially, may be because the correspondence between the 
receptor conformational state and CheA activity is more complicated, in connection with r\ 
and r_i discussed above. 



IV. SUMMARY AND DISCUSSIONS 

We suggest that statistical mechanics is useful and important in understanding receptor 
signalling and adaptation. We have made semi-quantitative comparisons between the theory 
and recent experiments to obtain estimations of parameter values. The thermal fluctuation 
in a cell is very strong, ksT w 4pN ■ nm ~ 0.025eV, comparable to the energy scale. So 
we simplify the formulation by using high temperature approximation. Then Eqs. (0) and 
( flOf) essentially contain all the information we need. 1 — f3u J characterizes the enhancement 
of signalling by receptor-receptor interaction. With this simplified formulation, we look at 
recent experimental results. Unlike a clean system usually studied in physics, for such a 
complex system, we do not expect the fitting to be quantitatively perfect. From the data 
on pre-stimulus motor rotation bias ||23|| , we obtain the pre-stimulus activity, as in Eq. (|i~7|), 
implying that there are approximately 70% receptor dimers are at the state corresponding the 
lower rate of CheA autophosphorylation. Although the data on response-stimulus relation 
are very limited, they are used to estimate that /3B/(1 — (3vJ) ~ 10.49, which compares 
the effect of ligand binding with that of cooperation. We study adaptation time for two 



different sets of data [£2|,£(j]. Assuming the delayed time in feedback to be 0.1s, it is found 
that the feedback strength compared with coupling, (3a /{I — (3vJ), is approximately 0.0068 
to 0.013, or 0.0015 to 0.0047, respectively. These numbers obtained from different data and 
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by using different methods are of the same order of magnitude, as a sign of the consistency 
of theory. Precise information on the feedback delay time can improve this determination. 
From the data on the relative CheA activity |P^| , we obtain Eq. (|27|), which gives the 
relation between the two levels of CheA activity corresponding to the two conformations of 
the receptor dimer. Combined with other results, it tells that the ratio between the two levels 
of CheA activity is A-i/A\ ~ 164.77, in good consistency with the available experimental 
information on this ratio. We note that the fitting is not perfect. This may be partly due to 
the simple nature of the minimum model and further simplified treatment, and partly due 
to insufficient experimental information. However, with a working framework proposed, we 
anticipate more experimental and theoretical discoveries stimulated by the current attempt. 
On the other hand, it would not be satisfactory to us to have a good fitting of the data by 
simply tuning parameters without a clear physical picture. 

We need improvement on the available experimental results, as well as new experimental 
information, to provide a basis for extension and refinement of the theory. For example, we 
need a significant broadening of the range of occupancy change in response-stimulus relation. 
We also need a clearer relation between adaptation time and occupancy change. The relation 
between CheA activity and the receptor conformational state and CheA activity, as well as 
the relative rate of the two rotation modes, is vital for going beyond the simple treatment 
here. More accurate results on A_\j A\ is also important. Independent determination of 
the dissociate constant is of fundamental importance. Most exciting experiments might be 
direct measurements of the conformational states V°, V 1 , and the coupling coefficient T^-, 
as well as the energy change or force induced by the ligand binding. A clarification on 
whether the conformational change is rotation or a vertical displacement is interesting. For 
the former, V° and V 1 are angles, while H, the effect of ligand binding, is a torque. For the 
latter, V° and V 1 are positions, while H is a force. The receptor-receptor interaction can 
be determined by measuring the relation of force or torque on one receptor dimer and the 
conformations of its neighbors. This would be a direct test of the conformation-dependent 
interaction. A determination of the geometry of the lattice is also interesting, from which 
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one can obtain the value of (3vJ, and consequently other parameter values. 

Our theory is entirely different from Hill model. An integer Hill coefficient is understood 
as the number of ligands bound to a receptor. A non- integer Hill coefficient, as often 
used, does not seem clear conceptually though could be tuned to fit the data. Nonetheless, 
from mean field point of view, the effect of receptor-receptor interaction could be viewed 
as effective additional ligand binding. Therefore from this perspective, the conclusion on 
limited cooperativity in [14| is consistent with strong thermal fluctuation in our theory. 



Here we specialize in chemotactic receptors, however, the theory may also apply to many 
other receptor systems. For example, state-dependent co-inhibition between transmitter- 
gated cation channels was observed p4|. Clustering of GAB A a receptors and the decrease 



of affinity was also studied p5|, which was also analyzed in terms of Hill model in a sim- 



ilar way to fl4| , thus it can also be explained by our theory as an indication of receptor- 
receptor interaction and thermal fluctuation. In many receptor systems, clustering, or called 
oligomerization, together with signalling, occurs as a response to stimulus. Theoretical in- 
vestigation on this situation is presented elsewhere. 

In finishing this paper, let us list some new experiments anticipated from the point of view 
of this theory. (1) More clarifications on conformational change induced by ligand binding, 
and determination of conformational change due to interaction with another receptor dimer. 
(2) Direct measurement of the forces generated by ligand binding and by conformational 
change of the neighboring receptor dimer. (3) Independent determination of dissociate 
constant using other methods. (4) Investigations on the responses corresponding to fixed 
conformational states, and hence r\ and r_x as discussed above. (5) Direct measurements 
on CheA and CheY activities. (6) More clarification on the relation between the receptor 
conformational state and CheA activity. (7) Increasing the range of occupancy change in 
response-stimulus relations, and more accurate determination of pre-stimulus occupancy 
and occupancy change. (8) More accurate determination of adaptation time as a function 
of the occupancy change. (9) Precise determination of the form of energy function. (10) 
Determination of the details of feedback due to the change of the methylation level, including 
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the delay time. 



17 



Figure captions: 

Fig. 1. An illustrative snapshot of the configuration of receptor dimers on a 50 x 50 
square lattice. An up triangle represents the conformation state Si — 1, a down triangle 
represents Si = —1, a filled triangle represents ligand binding, an empty triangle represents 
no ligand binding. 

Fig. 2. Response-stimulus relation 5m vs. Sc. The data points are transformed from 
those read from jl3| with a computer software. The range of receptor occupancy change is 
too small, so only qualitative comparison is possible. The straight line is the least square 
fitting Sm = 10.495c. 

Fig. 3. (a) Normal-normal plot of the relation between adaptation time t* and occupancy 
change Sc. The data points are adopted from |22|] by using a computer software, with the 
concentration transformed to occupancy, (b) Normal-log plot of the same data, showing 
that they can be fitted to a logarithmic relation. 

Fig. 4. Relation between adaptation time t* and occupancy change Sc. The data points 
are adopted from [2D|| by using a computer software, with the concentration transformed to 
occupancy. The straight line is the least square fitting t* = 156.3513 log 10 5c + 114.9912. 

Fig. 5. The relation between the relative CheA autophosphorylation rate R and ligand 
concentration L. The data points are adopted from fl4| by using a computer software. The 
theoretical curve is R = 1 — F L + Kd , with F = 0.95 and = 20/iM. 
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